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Abstract 

Nonresonant three-body decays of D and B mesons are studied. It is pointed out that if 
heavy meson chiral perturbation theory (HMChPT) is applied to the heavy-light strong and weak 
vertices and assumed to be valid over the whole kinematic region, then the predicted decay rates for 
nonresonant charmless 3-body B decays will be too large and especially B~ — * ir~K + K~ greatly 
exceeds the current experimental limit. This can be understood as chiral symmetry has been 
applied there twice beyond its region of validity. If HMChPT is applied only to the strong vertex 
and the weak transition is accounted for by the form factors, the dominant B* pole contribution 
to the tree-dominated direct three-body B decays will become small and the branching ratio will 
be of order 1(U 6 . The decay modes B~ — > (.fr~/i + /i~ )nr and 2?° — ► (If °/j + /j~)nr for h = ir,K 
are penguin dominated. We apply HMChPT in two different cases to study the direct 3-body D 
decays and compare the results with experiment. The preliminary FOCUS measurement of the 
direct decay Df — ► (vr + 7r + 7r _ )NR may provide the first indication of the importance of final-state 
interactions for the weak annihilation process in nonresonant D decays. Theoretical uncertainties 
are discussed. 
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I. INTRODUCTION 



The three-body decays of heavy mesons are in general dominated by intermediate (vector 
or scalar) resonances, namely, they proceed via quasi-two-body decays containing a resonance 
state and a pseudoscalar meson. The analysis of these decays using the Dalitz plot technique 
enables one to study the properties of various resonances. The nonresonant contribution is 
usually a small fraction of the total 3-body decay rate. Nevertheless, its study is important for 
several reasons. First, the interference between resonant and nonresonant decay amplitudes 
in B decays may provide information on the CP-violating phase angles |]-||. For example, 
the interference between B — > (7r + 7r~7r~)NR and B~ — > Xco 71 ^ could lead to a measurable 
CP asymmetry characterized by the phase angle 7 |J, while the Dalitz plot analysis of 
B — > pn — > 7T7T7r allows one to measure the angle a. Second, an inadequate extraction of the 
nonresonant contribution could yield incorrect measurements for the resonant channels [0. 
Third, some of nonresonant 3-body D decays have been measured. It is thus important to 
understand their underlying mechanisms. Experimentally, it is hard to measure the direct 3- 
body decays as the interference between nonresonant and quasi-two-body amplitudes makes 
it difficult to disentangle these two distinct contributions and extract the nonresonant one. 

The direct three-body decays of mesons in general receive two distinct contributions: 
one from the point-like weak transition and the other from the pole diagrams which involve 
four-point strong vertices. For D decays, attempts of applying the effective SU (4) x SU{4) 
chiral Lagrangian to describe the DP — > DP and PP — > PP scattering at energies ~ 
have been made by several authors [^-[T2| to calculate the nonresonant D decays, though 
in principle it is not justify to employ the SU(4) chiral symmetry. As shown in fTI|]l2]| , the 
predictions of the nonresonant decay rates in chiral perturbation theory are in general too 
small when compared with experiment. 

With the advent of heavy quark symmetry and its combination with chiral symmetry 
T3| |i~5| , the nonresonant D decays can be studied reliably at least in the kinematical region 



where the final pseuodscalar mesons are soft. Some of the direct 3-body D decays were 
studied based on this approach [fT~6|JT7[1 . 

Nonresonant charmless three-body B decays have been recently studied extensively based 
on heavy meson chiral perturbation theory (HMChPT). However, the predicted decay rates 
are unexpectedly large. For example, the branching ratio of B~ — > (7r + 7r~7r~)NR is predicted 
to be of order 10~ 5 in |I| and 0. Therefore, it has a decay rate larger than the two-body 
counterpart B — > tttt. However, it is found in || that the dominant B* pole contribution 
to the nonresonant B~ — > tt + 7t~tt~ accounts for a branching ratio of order only 1 x 10 -6 . 
Recently, Belle [18] and BaBar [19] have measured several charmless three-body B decays 
without making any assumptions on the intermediate resonance states |18[]. The predicted 
branching ratio of order 3 x 1CT 5 in [[| for B~ — > (K~ K + 7r~) NR already exceeds the upper 



limit 1.2 x 10 by Belle [18] and 7 x 10 by BaBar [[L9| for resonant and nonresonant 



contributions. Likewise, the predicted B{B — > it + it n )nr ~ 4 x 10 5 in is too large 
compared to the limit 1.5 x 10~ 5 set by BaBar. Therefore, it is important to reexamine and 



2 



clarify the existing calculations. 

The issue has to do with the applicability of HMChPT. In order to apply this approach, 
two of the final-state pseudoscalars have to be soft. The momentum of the soft pseudoscalar 
should be smaller than the chiral symmetry breaking scale A x ~ 830 MeV. For 3-body 
charmless B decays, the available phase space where chiral perturbation theory is applicable 
is only a small fraction of the whole Dalitz plot. Therefore, it is not justified to apply chiral 
and heavy quark symmetries to a certain kinematic region and then generalize it to the region 
beyond its validity. In order to have a reliable prediction for the total rate of direct 3-body 
decays, one should try to utilize chiral symmetry to a minimum. Therefore, we will apply 
HMChPT only to the strong vertex and use the form factors to describe the weak vertex. 
In contrast, for direct 3-body D decays, the allowed phase space region where HMChPT is 
applicable can be a dominant one for some decay modes. 

The paper is organized as follows. After introducing the effective Hamiltonian in Sec. II 
we proceed to discuss the difficulties with HMChPT when applying it to describe the 3-body 
nonresonant B decays in the whole Dalitz plot and its possible remedy. The full amplitude 
for the penguin-dominated B~ — > K~ir + ir~ is worked out as an example. The direct 3-body 
D decays are discussed in Sec. III. Discussions of theoretical uncertainties and conclusions 
are presented in Sec. IV. 



II. NONRESONANT THREE-BODY DECAYS OF B MESONS 

A. Hamiltonian 

The relevant effective AB = 1 weak Hamiltonian for charmless hadronic B decays is 



n eS (AB = i) = ^ Vub v: q \cMo u M + cMo u m] + v cb v; Q \ Cl ^)oi^) + c 2 (/i)o 2 » 



uq 

io >| 

VaV^CiitiO^) +h.c, (2.1) 

i=3 ) 



V-A J 
V-A 1 



where q — d,s, and 

°i = (ub) v _ A (qu) v _ A , 01 = {u a b j3 ) v _ A {qpu a ) 

0{ = (cb) v _ A (qc) v _ A , O c 2 = (c a b l3 ) v _ A (q p c a ), 

03(5) = (Q b )v-A £(9V)v-a(v-m)> °4(6) = (qah)v-A HW*)v-aW+a)i ( 2 - 2 ) 

q> q> 

3 3 

°7(9) = ^{q»v- A Y. e ^^ < l)v + AW-A)^ 08(io) = 2^ ab ^ v - A Y. e l'(ma)v + A(v-A)i 

q' q' 

with O^-Oq being the QCD penguin operators, O7-O10 the electroweak penguin operators 
and {qiq2) v±A = 5i7ai(1 ^ 15)0.2- The scale dependent Wilson coefficients calculated at next- 
to-leading order are renormalization scheme dependent. In the factorization approach the 
decay amplitude has the form 
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A(B -> MiM 2 M 3 ) oc 2 o< ( M i M 2M 3 1 5), (2.3) 

where the coefficients aj are renormalization scale and 75-scheme independent. In ensuing 
calculations we will employ the values of a, listed in |20||. For D decays we will use 



ai = 1.20, a 2 = -0.67. (2.4) 



B. Difficulties with heavy meson chiral perturbation theory for nonresonant B decays 

The nonresonant three-body B decays have been studied in two distinct methods, though 
both are based on heavy quark symmetry. One relies heavily on chiral perturbation theory to 
evaluate the 3-body matrix elements Pf3|,|21"| , whereas the use of chiral symmetry is restricted 



to the strong vertex for the other case [|TJ,^J| . The resulting decay rates can be different by 
one to two orders of magnitude. 

Let us first recapitulate the approach of heavy meson chiral perturbation theory P^-IB| 



and consider the decay mode B — ► (K K + n )nr as an illustration. Since this decay is tree 
dominated, we will focus on the dominant contribution from the four-quark operator 0\ 

A{B- - K~( Pl )K + (p 2 )n-(p 3 )) = ^VrtV^^iK-K+n-lO^B-). (2.5) 

Under the factorization approximation, 

{K-K+ir-\O x \B-) = (T-\(du) v _ A \0)(K-K + \{ub) v _ A \B-) 

+ (K-K+7T-\(du) v _ A \0)(0\(vh) v _ A \B-). (2.6) 

The second term on the right hand side corresponds to weak annihilation and it is expected 
to be helicity suppressed. As we shall see below, it indeed vanishes in the chiral limit. 
The three-body matrix element {K~ K + \(ub) v _ A \B~) has the general expression 



(K (p 1 )K + (p 2 )\(ub) v _ A \B (p B )} = ir(p B - Pi -paV + w+(p2 +P1V 

+ iu_(p 2 -P\)p + he^ af3 p B (p 2 + Pi) a (p 2 - Pi) 13 , (2.7) 

where r, u± and h are the unknown form factors. When pseudoscalar mesons are soft, the 
heavy-to-light current in the heavy quark limit can be expressed in terms of a heavy meson 
and light pseudoscalar mesons p4]JT3[] . The weak current = g a 7^(l — 15)Q, when written 
in terms of a heavy meson and light pseudoscalars, has the form 



K = %fHb f^ Tr[ 7 "(l - l5 )H h i\ a ] (2.8) 



to the lowest order in the light meson derivatives, where H a contains the pseudoscalar meson 

a/i' 



P a and the vector-meson field P*„ 
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FIG. 1. Point-like and pole diagrams responsible for the B — > K K + matrix element of the 
current ^7^(1 — 7s)fr, where the symbol • denotes an insertion of the current. 



Ha = V^^iP^ - Pal,) 



(2.9) 



where v is the velocity of the heavy meson and £ 2 is equal to the unitary matrix U which 
describes the Goldstone bosons. The general expression of the matrix U up to the fourth 
order in the meson matrix 6 is |23| 



U = 1 + 2if - 2^ - ia 3 ^ + 2(a 3 - 1)^ + • 

Jt Jtt Jit Jit 



(2.10) 



where a 3 indicates the nonlinear chiral realization and it has the well-known value | in the 
usual exponential expression for U, namely, U = exp(z20// 7r ). Here we do not specify the 
value of a 3 in order to demonstrate that the physical quantity is independent of the choice 
of chiral realization, i.e. the value of a 3 . The traceless meson matrix reads 



/ 2£L + JL 



7T 



71 



V2 + 



JL 



K+ 
K° 



V K- 



K 



(2.11) 



To compute the form factors r, u± and h, one needs to consider not only the point-like 
contact diagram, Fig. 1(a), but also various pole diagrams shown in Fig. 1. The heavy 



meson chiral Lagrangian given in |[T3Hl5|] is needed to compute the strong B*BP, B*B*P 
and BBPP vertices. The results for the form factors are [p2],2| 



g fB;m B . y /m B m B . 

fl t - m 2 B . 
g fB*m Bty /m B m Bt 
fl t - m 2 B , 



1 - 



(Pb - Pi) - Pi 



m 



+ 



B 



1 + 



(Pb - Pi) "Pi 



ml 
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r = £b_ _ h_ Pb ■ (P2-P1) 2gf B * j m B (pb-Pi) - Pi 

VI fl (Pb - Pi - P2) 2 - m% fl Y m B * t-m B% 

_ 4g 2 f B m B m B . Pi -p 2 - Pr (pb - Pi) P 2 - (gg - Pi)/m% ^ _ 

/I (Pb - Pi - P2) 2 -m% t- m| ; 

with £ = (pb — Pi) 2 = (P2 +P3) 2 - Note that the term f B /(2f%) comes from the point-like 
diagram, while the other terms in uj + and u;_ arise from the B* pole contributions in Fig. 1. 
The decay amplitude then reads 

A(B~ - K-( Pl )K + (p 2 )n-(p 3 )) NR = _^KdK>iy (2-13) 



x 



{2m 2 r + (m B — s — m 2 )uj + + (2t + s — m B — 2m\ — mg)u;_| , 



with s = (p B — P3) 2 = (pi +P2) 2 - It is clear that the contribution due to the form factor r is 
proportional to m 2 and hence negligible. For the strong coupling g, which will be introduced 
again below, we shall employ the value of g = 0.59 ± 0.01 ± 0.07 as extracted from the recent 
CLEO measurement of the D* + decay width [^J] . 

The decay rate of B~ — > K~K + it~ is then given by 



r (B - K-K+*-) = (^3^/^ \A?isM. (2.14) 

For a given s, the upper and lower bounds of t is fixed. If (|2.13|) is applicable to the whole 
kinematical region, then s min = (mi +m 2 ) 2 and s max = (m B — m 3 ) 2 , and the branching ratio 
of B~ — > K~ K + ix~ is found to be 

2.8 x 10 -5 from the contact term only, 
-> ^K + 7r~) NR = <( 6.7 x 10" 5 from the B* pole only, (2.15) 

1.7 x 10" 4 total. 



This is already above the upper limit of 7.5 x 10 5 set by CLEO PSl, and it greatly exceeds 



the experimental limit 1.2 x 10 reported recently by Belle [18j and 7 x 10 by BaBar 



T9| , recalling that both Belle and BaBar do not make any assumptions about intermediate 
resonances. In other words, the upper bound on the nonresonant B~ — > tt~K + K~ is pre- 
sumably much less than 1 x 10~ 5 after subtracting resonant contributions. Therefore, it is 
very likely that the branching ratio of direct B — > PPP decays is overestimated by one to 
two orders of magnitude in this approach. 

The dominant contributions to the direct B~ — ► K~K + tt~ come from the B* pole and 
the point-like weak transition term fs/f 2 - Since the chiral representation for the heavy-to- 
light current is valid only for low momentum pseudoscalars, the contact contribution from 
(7T~\(du)\0)(K + K~\(ub)\B~) and the weak B* to K transition in the B* pole diagrams are 
reliable only in the kinematic region where K + and K~ are soft. Therefore, the available 
phase space where chiral perturbation theory is applicable is very limited. It is claimed in 
PJ3|J21[| that if the usual HQET Feynman rules for the vertices near and outside the zero- 
recoil region but the complete propagators instead of the usual HQET propagator are used, 
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then the model is applicable to the whole Dalitz plot. However, as shown above, this will lead 
to too large decay rates in disagreement with experiment. Therefore, in order to estimate the 
nonresonant rates for the whole kinematic region, one should try to apply chiral symmetry 
to a minimum or some assumptions have to be made to extrapolate chiral symmetry results 
to the whole phase space. 

C. B* pole contribution 

As discussed before, the direct contact contribution to the matrix element 
(K + K~\(ub) v _ A \B~) as characterized by the fslfl term is valid only in the chiral limit, 
and hence we will not consider its contribution when computing the total decay rate. As for 
the B* pole contribution, we shall try to avoid the use of chiral symmetry when computing 
the B* to K weak transition; that is, we shall not use Eq. ( |2.8| ) to evaluate the matrix 
element of the B* — > P transition and we apply HMChPT only to the strong vertex and use 
form factors to describe the weak vertices. In this way, the soft meson limit is applied only 
once rather than twice. 

For the tree-dominated decay B~ — > K~K + -jt~ , the B* pole contribution isf] 

A B%vK 2 _ 2 BB*K' l Z - iD J 
¥B* B* 

The general expression for A U BB * K is 

e v A BBtK = (K-(q)B + (p B )\Bf(p B *)) = g BB * K (e ■ q). (2.17) 



In heavy quark and chiral limits, the strong coupling Qbb%k is determined to be |T3|-^5|1 



2g , 

9bb*k = -ry/m B m B *, (2.18) 

Jtt 



where g is a heavy-flavor independent strong coupling and its sign is positive ||13|| . It should 
be stressed that the relation ( |2.18| ) is valid only when the kaon is soft. Under the factorization 
approximation 

e,A B>K = ^K^«i(^(P3)|(^) v _ A |0)(^ + (p 2 )|(^) v _j5:°). (2.19) 

Heavy quark symmetry is then applied to relate the matrix element of B s — > K + to B s — > 
K + ffl: 



*The pole contribution from the scalar meson Bq and the effect of the decay width in the propagator 
have been considered in 0. We find these effects are small. 
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(K + (p K )\(ub) v _ A \B* s °(p B *)) = T x ie» va pe v VB*VK ~ T 2 Tr? Bt e il - T 3 (e ■ Pk)(pb* s + Pk)^ 

- T±(£ ■ p K )(p B * -p K )^ 
(K + (p K )\(ub) v „ A \B° s (p Bs )) = U(p Bs +p K ), + f-(pB 3 -p K )», (2.20) 

with being the polarization vector of B* s . The result is (see e.g. [|T 



T 3 



/+-/- 

m B 

/+-/- 



m 



2m 



B 



(/+ + f-)m B + (/+ ~ /-) 



Pb* ■ Pk 
m B 



B 



In terms of the form factors F^q defined by 



(2-21) 



2 2 

B a K/2\ , m B s ~ m K n \ TpB a K i 2 



(K + (p K )\(ub) v _M(PB)) = (PB+PK)»Ff° K (q 2 ) + 



q,[F BsK V) ~ F^ K \q 2 )} (2.22) 



with g M = (p B - Pk)^ we obtai 



am 



/+ = n 



B a K 



f- 



rn 



B TpB s K 



mi 



1 - 



7 B S K- 




Ft 



B a K ' 



and 



.G f 



t 

m B H m B - 

m B 



m B — t 



mi 



Hence, the B* pole contribution to B — ► K K + n is 



A(B~ - ^-(pi)^ + (p 2 )7r-(p 3 ))poie = %V ub V: dai F^ K {ml)-^ z 

\/ A i — m B » 



x 



X 



t m B — t 
m B H m B — 

m B mi. 



1 - 



F B ° K (ml) 
F*° K {ml) 



(2.23) 



(2.24) 



(2.25) 



' , f TO a 2 {t-ml + m\){m 2 B -t-m\) 



2m\, 



Using the Melikov-Stech model [27]] for the B s ^ K form factors, the branching ratio due 
to the B* pole is found to be of order 1.8 x 10~ 6 , which is consistent with the upper limit 



s 
1-5 



1.2 x 10" 5 set by Belle pfl and 7 x l(T b by BaBar [IS 



tit is most convenient to apply the interpolating field method for heavy mesons (see e.g. 
namely, \B*) = h^^q and \B) = h^i^q, to relate the B* — > P form factors to those of B — > P. 
The matrix element (ir + \(ub) v _ A \B°) is also evaluated in [|| using the relativistic potential model. 
However, only the form factor T 2 is calculated there. 
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In contrast, the matrix element of B — > in HMChPT has the form 



(K + {p K )\{ub) v _M{PBt)) 



■m B *e 



(2.26) 



Comparing this with Eqs. ( |2.20| ) and ( p.21[ ) it is clear that in the heavy quark and chiral 
{Pk 0) limits, only the form factor T 2 contributes with 



m B T 2 



f Bt f * 



fir fir 



— - in heavy quark and chiral limits, 



(2.27) 



where use of Eq. ( |2.23| ) has been made. However, beyond the chiral limit, all T 2 , T 3 and T4 
contribute and 



m B T 2 = F^ K {ml\ 



1 + 



t 



m: 



+ m\ 



2m| 



t + mt 



m 



^B S K 



K 



2m 2 B 



2ml 



1 - 



(2.28) 



in the heavy quark limit. Since F^ aK (0) = 0.31 in the MS form-factor model ||27|| , it is 
evident that the form factor T 2 inferred from Eq. ( [2.28P is much smaller than that implied 
by Eq. ( |2.27| ), namely, T 2 = f B J ' fn = 1-6 for f Bs = 190 MeV. This explains why the 
prediction based on HMChPT is too large by one to two orders of magnitude compared to 
the B* pole contribution which relies on chiral symmetry only at the strong vertex. 

The previous estimate of B~ — > (7r + 7r~7r~)NR by Deshpande et al. based on the B* 
pole contribution gives a branching ratio of order 2 x 10 -5 for Ff 71 ^) = 0.333 and g = 0.60 
(case 1 in []]]). This is larger than our result 3.0 x 10~ 6 (see Table I below) by one order of 
magnitude. It can be traced back to the square bracket term in Eq. ( |2.24| ) for the analogous 



£ fA B * 7r7T 



term where Deshpande et al. obtained 



3 

-m B + 
2 2m B 



t m B m B — t 



mi 




(2.29) 



to be compared with 



m B + 



m. 



m B 



m B - 



1 - 



mi 




(2.30) 



in our case. Numerically, the decay rate obtained by Deshpande et al. is larger than ours 
by a factor of 3 when the same B — > tt form factors are employed. Note that the B* pole 
contribution to B~ — > 7r + 7r~7r~ is found to be 1.8 x 10 -6 (for g = 0.6) in and 2.7 x 10~ 6 
in ||. Therefore, our result is consistent with them. 

D. Full contributions 



In the previous subsections we have only considered the dominant contribution to the 
tree-dominated B decay from the operator 0\. In the following we discuss the full amplitude 
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for the direct 3-body B decay and choose the penguin-dominated decay B — > n n + K as 
an example. The factorizable amplitude reads 



A(B~ - n-( Pl )n + (p 2 )K-(p 3 )) = %(v^C 



ai (K-\(su) v _ A \0)(n + n-\(ub) v _ A \B-) 



+ (n-n+K- 1 (su) v _ A \0) (0| (ub) v _ A \ B~) 
+ a 2 (n-n + \(uu) v _ A \0)(K-\(sb) v _ A \B-) 



(2.31) 



+ -(a 7 + a 9 )(n n + \(e u uu + e d dd) v _ A \0)(K \(sb) v _ A \B ) 



- v tb v* 



a 4 (rt-7i + K-\0 4 \B-) + a^ir"^ R-p^B') 



+ (4 -> 10) + (6^8) 



Under the factorization approximation, the matrix element of O4 is 

(Tr-TT+iT-lCMS-) = (^|(^) v _jO)(7T-7r + |(n6) v _J J B-) 
+ {7r + K-\(sd) v _ A \0){n-\(db) v _ A \B-) 

+ (Tr-Tr+jr-Kfuj^ioxoK^^is-). 



(2.32) 



In Eq. ( |2.31|) the two-body matrix element (n + K \(sd) v _ A \0) has the form 



(n + (p 2 )K-(p 3 )\(sd) v „ A \0) = (n + (p 2 )\(sd) v _ A \K + (-p 3 )) = (p 3 - P2) ^ (t) 



+ 



rn 



K 



m; 



-(P3 + P2), -F?*{t)+Ff*{t) 



(2.33) 



where we have taken into account the sign flip arising from interchanging the operators 
s <-> d. The other two-body matrix element (tt + 7t~\(uu) v _ a \0) can be related to the pion 
matrix element of the electromagnetic current 



(2.34) 



with q 2 = {p' — p) 2 and J^ m = ^wy^u — \d^^d + ■ ■ ■. The electromagnetic form factor F n7T is 



normalized to unity at q 2 = 0. Applying the isospin relations yields 

(7r + (p)\u^u\7r + (p')) = (7r-(p)\d llx d\7r-(p')) = {p + p'),F™{q 2 ). 



(2.35) 



As for the three-body matrix element (tt~tt + 'K~\(su) v _ a \0) , one may argue that it van- 
ishes in the chiral limit owing to the helicity suppression. To see this is indeed the case, we 
first assume that the kaon and pions are soft. The weak current can be expressed in terms 
of the chiral representation derived from the chiral Lagrangian 



x2 x2 

£ = J -^Tr(d^Ud^) + J -^-Ti(MU ] + U ] M). 
8 ' 8 



(2.36) 
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The weak current = %7 M (1 — 75)A a (^ has the chiral representation (see e.g. 



J" = -Ijl Tr{U^X a d^U - d^X a U) = Tr^A^t/). 
It is straightforward to show that J M = ^7^(1 — 75)^ has the expression 

J ] ; = -4- i^r 9 ^ + M - ^0} + ^(4 - a 3 )^ + 



(2.37) 



31 



P 

J 7T 



P 



f 3 



Note that the sign convention of J® or J M is chosen in such a way that (0| J^\P{p)) 
We are ready to evaluate the point-like 3-body matrix element 



(n-( Pl )7t+(p 2 )K-(p 3 )\(su) v „ A \0) c 



1 



a 3 



(Pi + P2 +P3) M - 2p : 



2/t 



(2.38) 
(2.39) 



which is chiral-realization dependent. This realization dependence should be compensated 
by the pole contribution, namely, the B~ to K~ weak transition followed by the strong 
interaction K~ — > K^-k + -k^ . The strong vertex followed from the chiral Lagrangian ( [2.36|) 
has the form 



S 



— (p -m 3 ) + — p-p 2 , 

Z /tt Jtt 



with p = pi + p 2 + P3. Hence, 

(n-( Pl )ir + (p 2 )K-(p 3 )\(su) v _ A \0) = K-\(su) v _ A \0) contiiCt + S 

P-P2 



p 2 — m 2 K 



2i 

In 



P2n 



> o P^ 

p z — mr K , 



(2.40) 



(K-{p)\(su) v _ A \0) 



:2AV 



Evidently, the 03 terms are cancelled as it should be. It is worth stressing again that the 
above matrix element is valid only for low-momentum pseudoscalars. It is easily seen that 
in the chiral limit 



(n-n+K-\(su) v _ A \0)(0\(ud) v _ A \B-) = 0. 



(2.42) 



Physically, the helicity suppression is perfect when light final-state pseudoscalar mesons are 
massless. Although Eq. (|2.42|) is derived for soft Goldstone bosons, it should hold even for 
the energetic kaon and pions as the helicity suppression is expected to be more effective. 
The factorizable contributions due to the penguin operator 6 is 



(n-n + K-\0 6 \B-) = — 2{<AT _ |s(l + 75HO) (7^+^(1 - ^)b\B~) 
+ (7r + ^-|s(l + 7 5 )rf|0)(7r-|J(l-7 5 )6|fi-) 
+ (7r-7r + K-|s(l + 7 5 )w|0)(0|M(l-7 5 )6|fi-)}. 



(2.43) 



Applying equations of motion we obtain 
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2 

(ir-| S ll + 75)n|0)(7r-7r + |S(l-75)&| J B-) = ^(^-|(^) v+4 |0)(7r-7r + |(^) v+A |i?-} 

m b m s 

= ^{K-\{su) v ^)^-^ub) V - A \B-), (2.44) 

and 

(7r + (p 2 )K-(p 3 )\s(l + 75 )d|0)(7r-(p 1 )|d(l - T5 )&|^-> 

(P2 +PsY , +, w.-/ M^. „.|n\ vBnt 



-(n+(p 2 )K-(p 3 )\s^u\0) * Ft(t) 



2 

^F^^F^t). (2.45) 



To evaluate the three-body matrix element (tt~tt + K~\s(1 + 7s)w|0), we will first consider 
the case that the kaon and pions are soft and then assign a form factor to account for their 
momentum dependence. At low energies, it is known that the light-to-light current can be 
expressed in terms of light pseudoscalars (see e.g. |2c 



g,(l-7 5 ) gi = ^% ( 2 - 46 ) 



to the lowest order in the light meson derivatives, where 



ml+ m 2 K+ m\ — m\ 



m u + m d m u + m s m s - m d 



(2.47) 



characterizes the quark-order parameter (qq) which spontaneously breaks the chiral symme- 
try. It is easily seen that the point-like contact term yields 

(7^ + ^-1575^0) contact = ^J"- ( 2 -48) 

* J IT 

As before, this chiral-realization dependence should be compensated by the pole contribution, 
namely, the weak transition of B~ to K~ followed by the strong scattering K~ — ► K^-k + -k^ . 
Hence, 

(7r-(pi)7r + (p 2 )^(P3)|s75M|0) = (7r~7r + iT~ |s7 5 w|0) C ontact + S— r (K~(p)\s'y 5 u\0} 



p 2 — m\ 



(249) 



Therefore, the 03 terms are cancelled. Note that, contrary to the {V — A)(V — A) case where 
the weak annihilation vanishes in the chiral limit, the penguin-induced weak annihilation 
does not diminish in the same limit. This is so because the helicity suppression works for 
the (V - A)(V - A) interaction but not for the (S - P)(S + P) one. 
Putting everything together leads to 
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(it-it + K-\0 6 \B-) 



—^—(K~ | (SU) V _ A \0) (7T-7T+ | (ub) v _ A \B-) + 

m b m s 



i 11 i 
m K — m^ m B — 



m. 



x 



Ft{t)Ft{t) 



B~ 



K 



P 

J IX 



1 



2pi ■ P3 



[m 



m B — m 



B) 



K , 



m b 
(2.50) 



where the form factor F K7T ' K is needed to accommodate the fact that the final-state pseu- 
doscalars are energetic rather than soft. The full amplitude finally reads 



A(B' 



71 7T + K 



)NR 



Gi 



(0 



V ub V*ai - VtbVZ a 4 + a w - 2(a 6 + a 8 ) 



in 



K 



m b m s 



x (K-\(7su) v _ A \0)(n-n + \(ub) v _ A \B-) 



+ 



V u bV* s a 2 - V tb V t * s -(a 7 + a 9 )j F» K (s)F™(s)(t - u) 

7T f 4-\ ( m B ~ m l)( m K ~ m l) 



- V tb V* (a 4 - -a 10 ) [Ff<(t)Ff*(t) 



+ F^(t)F^(t)(m B + 2m 2 r + m 2 K -2s-t 



(m 2 B - ml){m 2 K - ml) 



- (2a 6 - a 8 ) 



m b 

F»«(t)Ff*(t) - 



m. 



f 2 

■J 7T 



1 - 



2P1 • P3 \ nft 



— m 



Bj 



K , 



(2.51) 



where u = (p B — Pi) 2 ■ As noted in passing, we should only consider the pole contribution to 
the 3-body matrix element {it~it + \{ub) v _ A \B~) so that 



(ir-(p3)|(^) v „JO)(7r^(p 1 )7r + (P2)|(^) v ._Ji?-} 
*i { m K)~p ; 5 — m B H ™>b- 



pole 



f n t - m 2 B * 



m 



B 



KJ 



m B 



m 



K 



m 



K) 



X 



' , f ™2 2 (t-m| + m|)(m|-t-mj) 
s -f i m B m 2 -\ ^— y- 



(2.52) 



The decay amplitudes for other decays B~ — > it~ (K~)h + h~ and B K h + h~ have the 
similar expressions as Eq. ( |2.51|) except for — > 7r + 7r~7r~ and 5~ — > K + K~ K~ where one 
also needs to add the contributions from the interchange s t and put a factor of 1/2 in 
the decay rate to account for the identical particle effect. 



E. Results and discussions 



Before proceeding to the numerical results, it is useful to express the direct 3-body 
decays of the heavy mesons in terms of some quark-graph amplitudes [fTT|j29|] : T\ and T 2) the 
color-allowed external W-emission tree diagrams; C\ and C 2 > the color-suppressed internal W- 
emission diagrams; S, the VT-exchange diagram; A, the VT-annihilation diagram; V\ and V 2 , 
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the penguin diagrams, and V a , the penguin-induced annihilation diagram. The quark-graph 
amplitudes of various 3-body B decays B — > nh + h~ and B — > Kh + h~ are summarized in 
Table I. As mentioned in [IT], the use of the quark-diagram amplitudes for three-body decays 



are in general momentum dependent. This means that unless its momentum dependence is 
known, the quark-diagram amplitudes of direct 3-body decays cannot be extracted from 
experiment without making further assumptions. Moreover, the momentum dependence of 
each quark-diagram amplitude varies from channel to channel. 

To consider the nonresonant contribution arising from the pion and kaon electromagnetic 
form factors F nn and F , we follow |TJ to use the parametrization 

^nonres(9 ) = i o /™2 i ir< TZT ' (2.53) 

and employ r* = 200 MeV, and m* = 600 MeV for the pion and 700 MeV for the kaon. The 
momentum dependence of the weak form factor F (q 2 ) is parametrized as 

F K *tf) = i 2 f A ? T A ' (2 ' 54) 

1 - q 2 /A x + iT*A x 



where A x w 830 MeV is the chiral-symmetry breaking scale |p3 . Likewise, the form factor 
pK-K-K a pp ear j n g j n Eq. (|2.50|) is assumed to be 



^'V) = J-^- (2.55) 

The predicted branching ratios for direct charmless 3-body B decays are shown in Table 
I. The decays B~ — ► 7i~h + h~ are tree dominated and their main contributions come from the 
B* pole. In contrast, the decays B~ — ► (K~h + h~)-NR and 5° — > (K°h + h~)NR for h = n,K 
are penguin dominated. When h = ir, the main contribution comes from the 2-body matrix 
elements of scalar densities, namely, the second term on the right hand side of Eq. ( 2"^ 



while the contribution from the three-body and one-body matrix elements of pseudoscalar 
densities [the first term of Eq. (|2.43| )1 characterized by the term 2a§m 2 K /{morris) in Eq. 



( p.51| ) is largely compensated by the 04 term. 

Direct three-body charmless B ± decays have been searched for by CLEO [25] with 



limits summarized in Table III The decays B — > tx K + K , K K + K and B° 



K 7T + 7r _ , K K + K~ were measured recently by Belle p^8|j30|] and BaBar [TJJ but without any 
assumptions on the intermediate states. It is interesting to note that the limits 1.2 x 10~ 5 set 
by Belle and 7 x 10~ 6 by BaBar for tt~K + K~ (resonant and nonresonant) is improved over 
the previous CLEO limit 7.5 x 10~ 5 for the nonresonant one. Needless to say, it is important 
to measure the nonresonant decay rates by B factories and compare them with theory. 

In the estimation of direct 3-body decay rates we have applied the B*BP strong coupling 
given by Eq. ( 2.18|) and the B* — ► P weak transition beyond their validity. Needless to 



say, this will cause some major theoretical uncertainties in the calculations because the 
strong B*BP coupling is derived under heavy quark and chiral symmetries and hence the 
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FIG. 2. Quark diagrams for the three-body decays of heavy mesons, where Q denotes a heavy 
quark. 

momentum of the soft pseudoscalar should be less than A x . For the energetic pseudoscalar, 
the intermediate B* state is far from its mass shell. It is assumed in [[]]] that the off-shellness 
of the B* pole is accounted for by replacing the term ^rug* in Eq. ( |2.18| ) by (p 2 B*) 1 ^ and it 
is found that the branching ratios are reduced by (30 ~ 40)% for B~ — > tt~K + K~ , 7r + 7r~7r~ 
as shown in Table |]| while B~ — > K~h + h~ for h = n,K remain essentially unaffected. Using 
the measured branching ratios (55.6 ±5.8 ±7.7) x 10~ 6 and (35.3 ±3.7 ±4.5) x 10~ 6 by Belle 
0, (59.2 ±4.7 ±4.9) x 10" 6 and (34.7 ± 2.0 ± 1.8) x 10~ 6 by BaBar @ for B~ -> K~~tt + 7t~ 
and B~ — > K~K + K~ , respectively, in conjunction with the calculated results for direct 
3-body decays, the corresponding fractions of nonresonant components are found to be 4% 
and 3%, respectively. 



III. NONRESONANT THREE-BODY DECAYS OF D MESONS 

For nonresonant three-body D decays, the applicability of HMChPT should be in a better 
position than the B meson case. In Table II the maximum momentum p of any of the decay 
products in the D rest frame is listed. As stressed in [J03], D — ► KKK are the decay modes 
where HMChPT can be reliably applied since p there is of order 545 MeV which is below the 
chiral symmetry breaking scale. For other Ktttt and KKn modes, the regime of the phase 
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TABLE I. Quark-diagram amplitudes and branching ratios for nonresonant 3-body charmless 
B decays. The prediction Bl heoi is made for g B B*K(n) = 2 9 1 fir x (m B m B * ) 1/2 while the # t 2 heor 



accounts the off-shellness of the B* s ^ by letting g BB * ^k(tt) = tyffit x y rn B^jP B * j ■ Experimental 
limits are taken from [31]. 



Decay mode 


Quark-diagram amplitude 


B 1 

thfior 




B 2 

thoor 




^expt | 






K.fcV^VaCTi + Ci + -4) + V tb V t * d V2(,Vi +v 2 +v a ) 


3.0 x 10" 


-6 


1.7 x 10" 


6 


< 4.1 x 10- 


5 


-> TT-K+K- 


VubV^Tj. + Ci + A) + VttV;* (Pi + P 2 + Pa) 


1.8 x 1(T 


-6 


1.3 x 10 _ 


6 


< 7.5 x 10" 


-5 


—> K~n + n~ 


V„6V^(Ti + Ci + A) + V t(> V t * (Pi + P 2 + Pa) 


2.4 x 10" 


-6 


2.3 x 10" 


6 


< 2.8 x 10" 


5 


-* K-K+K- 


VvbVZsVZp! +d+A) + V th V t %V2(Pi +V 2 +Va) 


9.1 x 10- 


-7 


8.5 x 10" 


7 


< 3.8 x 10" 


5 


B — > K 7r T 7r 


V ub V* B Ci + V tb V t %(Ti +V 2 +V a ) 


2.1 x 10- 


-6 


2.1 x 10" 


6 










V ub V* s (Tl + Cl) + VttV'iV! +V 2 + Pa) 


1.2 x 10" 


-6 


1.2 x 10- 


6 









space where HMChPT is applicable is not necessarily small. 

The calculations for nonresonant three-body decays of the charmed mesons proceed in 
the same way as the B meson case and they are performed in the framework of HMChPT 
for two different cases: (i) HMChPT is applied to both strong and weak vertices, and (ii) it 
is applied only to the strong vertex and the weak transition is accounted for by form factors. 
These two different cases are denoted by B a and B b , respectively, in Table II. Here we would 
like to point out some interesting physics. First, consider the decay D° — > 7^°7r + 7r~ '. In 
HMChPT its amplitude is given by 

A{D o _» 7 r-(p 1 )lt ) (p 2 )-K + (p 3 )) = -^V cs V: d {a x A x + a 2 A 2 ), (3.1) 



with 



a i/tt f 2 / 2 2 \ / 2 2 2 \ 1 

Ai = — |2m 3 r + (m D — s — m 3 )to + + [2t + s — m D — 2m 2 — m 3 )uj_ j 



2 

A 2 = ^2m\r + {m 2 D — u — m\)uj + + (2t + u — m 2 D — 2m\ — m 2 )u>-} , (3.2) 

where the form factors r, uj + and cu_ have similar expressions as Eq. ( |2.12| ). Since d\ and a 2 
in D decays are opposite in signs [see Eq. ( |2.4j) 1, it follows that the decay rate is suppressed 
owing to the destructive interference, see Table II. 

However, when HMChPT is applied only to the strong vertex, the main contribution to 
7r + 7r comes from the D* + pole, namely, the strong process D° — > n D* + followed 
by the weak transition D* + — > K tt + . Since it is known that the interference in D + — > K tt + 
is destructive, naively it is expected that the same destructive interference occurs in the 
nonresonant D° — > T^tt+tt^ decay. However, this is not the case. The D* pole amplitude is 

FD* ~ m D* 

Now under factorization 
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e,A^ K = ^V cs V: d { a 1 (7r + (p 3 )IM v „jO)(K (p 2 )|( S -c) F _ A | J D* + (p D 0) 

+ a 2 



(^(p 2 )\(M) v _M(^(P3)\(nc) v _ A \D* + (p D ,)}}. (3.4) 



Applying heavy quark symmetry one can relate the form factors in (K°\(sc) v _ A \D* + ) to 
those in {K°\(sc) v _ A \D + }: 



{W(p K )\(sc) v _ A \D + (p D )) = f^(q 2 )(p D +p K ), + f^(q 2 )(pD-p K )». 
We obtain 



.Gt 



t 



e,A^ nK = - i ^V cs V: d (e-p 3 ) «i U (/+ + f-) DKm D + (/+ — f-) DK 

- a 2 f K [(/+ + / )" :: m n + (/+ - / f ' 



m D 



(3.5) 



(3.6) 



It is interesting to note that although the interference is destructive in D* + — > K 7r + , it 
becomes constructive in the process D° — > tt~D* + — > 7r _ 7r + i^ . We see from Table II that 
B b is indeed much larger than B a for .D — > ~K°7t + tt~ . 



TABLE II. Quark-diagram amplitudes and branching ratios (in percent) for nonresonant 
3-body D decays, where p (in units of MeV) is the largest momentum any of the products can 
have in the D rest frame. Heavy meson chiral perturbation theory is applied to both heavy-light 
strong and weak vertices for the theoretical prediction B a , while it is applied only to the strong 
vertex for B h . Form factors for D — > it and D — > K transitions are taken from [25] and experi- 
mental results from [31]. For the recent measurements of the nonresonant decays D + —* K~n + ii + , 
D° -> K°K+K- and D+ -> tt+vt+tt", see the text. 



Decay mode 



P 



Quark-diagram amplitude 



na pb to 

"theor "theor 



expt 





842 




0.03 


0.17 


see text 


-► K-ir+ir 


844 


K^ts^i+Ci) 


0.61 


0.28 


1 n5 +0 - 51 


K°K+K- 


544 


V^(T 2 + C 2 + £) 


0.16 


0.01 


0.55 ±0.09 


D+ -> K°ir+ir° 


845 


^K*^(r 1+ co 


1.5 


0.7 


1.3 ±1.1 


— > K~ir + ir + 


845 


K^^Tx+Ci) 


6.5 


1.6 


8.6 ±0.8 


— > 7T + 7T + 7r~ 


908 VLdV£>/2(Ti + Ci + .4 + Pi) + K fl V£x/2(7>i) 


0.50 


0.067 


0.024 ±0.021 


K-K+-K+ 


744 


V^t&M + Vi) + F US F* (Ti + d + £ ) 


0.48 


0.004 


0.45 ±0.09 


D+ -» K-K + it + 


805 




1.0 


0.69 


0.9 ±0.4 


—* 7r + 7r + 7r _ 


959 


V ud V* s V2(A) 






0.005 ± 0.022 
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The nonresonant decay D° — > (K K + K~)^ deserves a special attention for two reasons. 
First, it is the only Cabibbo-allowed direct 3-body mode which receives contributions from 
the external jy-emission diagram T 2 (see Fig. 2). Second, as noted in passing, HMChPT is 
presumably most reliable for this mode. Its factorizable amplitude has the form 



A(D° -> K-( Pl )K + (p 2 )K°(p 3 )) 



nr = ^V ud V;L(K + K\ud) v ^ A \0)(K~\(sc) 



\D°) 



+ a 2 (K°\(sd) v _ A \0)(K-K' 



[uc 



\D°) 



+ a 2 (K~K + K°\(sd) v _ A \0) (0\(uc) v _ A \D°) 



(3.7) 



where the three terms on the right hand side correspond to the quark diagrams T 2 , C 2 and £, 
respectively. Proceeding as before and neglecting the H^-exchange contribution in the chiral 
limit, we obtain 



A(D° - K-( Pl )K + (p 2 )K°(p 3 )) NR = ^gy^jojAi + a 2 A' 2 } 



where 



U \ t-m 2 m 2 



u 



A' 2 = \2m\r + {m 2 D — s — m\)uj + + (2t + s — m 2 D — 2m\ — m\)ujJ^ , 
when HMChPT is applied to both strong and weak vertices, or 



(3. 



(3.9) 



A! 1 = F? K {t)F KK {t){s-u) 



A' 2 = F( 



D S K,2 s9y/m D m D * 



(m 



K) 



m 2 Dt 



m D + 



t 



in 



x 



s + t — m D — m\ + 



m D 

(t — m\ + rn|)(m 



m D - 



D 



1 - 



m 



K 



D a K i 



K) 



ml 



2m% : 



(3.10) 



when HMChPT is applied only to the strong vertex. Again, the form factors r, u + and uj- 
in Eq. ( |3.9| ) have the similar expressions as Eq. ( |2.12j ). 

It is clear from Table II that the predicted branching ratio B a of 0.16% for D° — > 
(K° K + K~)nr works much better than B b , though the former is still too small compared to 
the experimental value (0.55 ± 0.09)% |JT|. This decay was also considered by Zhang JIB] 
within the same framework of HMChPT, but his result 2.3 x 10 -4 for the branching ratio, 
which is similar to the prediction 2 x 10 -4 based on chiral perturbation theory ||11|| , is smaller 
than ours by one order of magnitude. 

Some simple relations among different modes follow from the quark diagram approach. 
For example, neglecting the weak annihilation and penguin contributions and the phase 
space difference among different modes, it is expected that 
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B{D + 


— > 7T + 7T + 7r 
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B(D+ - 


— > K~7T^~7\ — 


) TVTT3 
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B(D+ - 


-> K-K+ti^ 
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B(D+ - 


-> 7T + 7r+7T _ 
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2' 




B(D+ - 


-> K~K + TT+ 


)nr 


1 r(D+) 


B(D+ - 


-> iT~7r+7r + 


)nr 


2r(L>+) 



(3.11) 



The above anticipation can be checked against the experimental results. It is easily seen 
that the measured D + — > (7T7r7r)NR is too small compared to the theoretical prediction. 
For example, the observation that (7r + K + K~)^ 3> (7r + 7r + 7r _ )NR m D + decays is rather 
unexpected. 

We see from Table II that the predictions for case (i) denoted by B a are generally larger 
than case (ii) denoted by B b except for the decay D° — > K°7r + 7r~. Contrary to the B meson 
case where the predicted rates in these two different methods can differ by one to two orders 
of magnitude, B a and B h in some of D decays differ only by a factor of 2. It is also evident 
that in general B b s give a better agreement with experiment for many of the direct 3-body 
D decays , whereas B a works better for D° — ► K°K + K~ and D + — > K~K + ir + , though the 
prediction of the former mode by HMChPT is still too small compared to experiment. As 
noted in the Introduction, the early predictions based on SU(4) chiral perturbation theory 
are in general too small when compared with experiment flTT| , [T12|j . 

There exist several new measurements of direct 3-body D decays in past few years: 
D° -> K-n+ir , K°ti + tt-, K°K + K-, D + -> tt+tt+tt-, R--n + Ti + and D+ -> tt+tt+tt-. The 
nonresonant branching ratio for the first mode is found to be (1.0 ± 0.1 ± OT^oii) x 10~ 2 by 
CLEO |32]| . Previous experiments |^3[ indicate that the decay D + — > K^-k + -k + is strongly 
dominated by the nonresonant term with (95 ±7)% [[H]]. However, a recent Dalitz plot 



analysis by E791 J34J reveals that a best fit to the data is obtained if the presence of an 
additional scalar resonance k is included. As a consequence, the nonresonant decay fraction 
drops from 95% to (13±6)%, whereas K7r + accounts for (48±12)% of the total rate. Therefore, 
the branching ratio of the direct decay D + — > K^-k + -k + is dropped from (8.6 ± 0.8)% to 
(1.2 ± 0.6)%. Likewise, It was found by the E687 experiment that the decay D + — > 7r + 7r + 7r~ 
is dominated by the nonresonant contribution with (60 ±11)% [35]. Again, the new Dalitz 



plot analysis by E791 |$6|] points out that half of the decays is accounted for by the scalar 
resonance a, whereas the nonresonant fraction is only (7.8 ± 6.0 ± 2.7)%. Consequently, 
B(D + — > 7r + 7r + 7r~)NR drops to (0.024 ± 0.021)%. Very recently BaBar has reported the 
preliminary result of the Dalitz plot analysis of D° — > T(°K + K~ [37]. Its nonresonant 
fraction is estimated to be (0.4 ± 0.3 ± 0.8)% and hence is negligible. 

As for the direct decay D° — > K 7i + tt^ , the 2000 edition of Particle Data Group (PDG) 



38| quotes a value of (1.47 ± 0.24)% for its branching ratio. However, it is no longer cited 
in the 2002 PDG [RTf] as no evidence for a nonresonant component is seen according to 
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the most detailed analyses performed in |]39] . This is also confirmed by a very recent CLEO 
measurement of this decay mode which gives (0.9 ±0.4^3^2)% f° r the nonresonant fraction 



The Cabibbo-suppressed decay — > (7r + 7r + 7r ) NR proceeds only through the 
lU-annihilation diagram. The early E691 measurement gives R = B{Df — ► 
7i + n + 7i-) NR /B{D+ -> 07T+) = 0.29 ±0.09 ±0.03 gl). However, it was found to be negligible 
by E791 |H and its branching ratio is quoted to be (5 ± 22) x 10~ 5 by 2002 PDG (see Table 
II). Recently, FOCUS has reported the preliminary result: the nonresonant fraction is mea- 
sured to be (25.5±4.6)% @. This corresponds to B(D+ -> 7r+7r+7r-) NR = (2.6±0.9) x lO" 3 . 
Although the short-distance VU-annihilation vanishes in the chiral limit, the long-distance 
one can be induced from final-state rescattering (see e.g. K|). Therefore, the observation 
of direct — > 7t + 7t + tt~ implies the importance of final-state interactions for nonresonant 
decays. 

IV. CONCLUSIONS 

We have presented a systematical study of nonresonant three-body decays of D and B 
mesons. We first draw some conclusions from our analysis and then proceed to discuss the 
sources of theoretical uncertainties during the course of calculation. 

1. It is pointed out that if heavy meson chiral perturbation theory (HMChPT) is applied 
to the heavy-light strong and weak vertices and assumed to be valid over the whole 
kinematic region, then the predicted decay rates for nonresonant 3-body B decays 
will be too large and especially B~ — ► tt~K + K~ exceeds substantially the current 
experimental limit. This can be understood because chiral symmetry has been applied 
twice beyond its region of validity. 

2. If HMChPT is applied only to the strong vertex and the weak transition is accounted 
for by the form factors, the dominant B* pole contribution to the tree-dominated 
direct three-body B decays will become small and the branching ratio will be of order 
10 -6 . The decay modes B~ — > (K^h + h~)^ R and B° — > (K°h + h~)i<sR for h = ir, K are 
penguin dominated. 

3. We have considered the use of HMChPT in two different cases to study the direct 3- 
body D decays. We found that when HMChPT is applied only to the strong vertex, the 
predictions in general give a better agreement with experiment except for the decays 
D° -> K-tt + tt°,K°K + K~ and D+ -> K~K+n + where a full use of HMChPT to the 
weak vertices gives a better description. The D* + pole contribution to D° — > K n + n^ 



*For previous theoretical estimates, see []17| and 
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proceeds through external and internal PF-emission diagrams with constructive in- 
terference. The experimental observation that (tt + K + K~)nr 3> (7r + 7T + 7r~)NR in D + 
decays is largely unanticipated. 

It is useful to summarize the theoretical uncertainties encountered in the present paper, 
though most of them have been discussed before: 

• For B* (and also D*) pole contributions, the intermediate state B* is off its mass shell 
when the pseudoscalar meson coupled to B* and B is no longer soft. This will affect the 
B*BP strong coupling. To estimate the off-shell effect of B*, we replace its mass mg. 
by \Jp 2 B , and find that the branching ratios for B~ — > n~K + K~ , 7r + 7r~7r~ are reduced 
by (30 ~ 40)%, while B~ — > K~tt + tt~, K~K + K~ remain essentially unaffected. 

• We have parametrized the q 2 dependence of the form factors F™ nres , F^ Tes , F^ Tes and 
pK-K-K - m f orm Q f gq S> (|2,53|) , (|2.54j) and (p.55[) . However, part of scalar resonance 



effects is included in the parametrization of the form factors. In the B decays, the 
major uncertainty of the calculated amplitudes comes from the chiral enhanced term 
~ F^ onres (2aQ — as) x ms(m|- —m%)/m s . We may overestimate the penguin-dominant 
nonresonant branching ratios if there exist scalar resonances, e.g. k. Although in some 
channels the a resonance is included in F^ me8 , its effect is suppressed by the Cabibbo 
angle and by the fact that it decouples to the vector current in the SU(2) symmetry 
limit. 

The point-like contact contribution to the three-body matrix element beyond the chiral 
limit, e.g. (PiP2\(qb) v _ A \B) contaiCt , is unknown but it becomes even smaller when P 1 
or P2 is not soft owing to the smaller wave function overlap among P\,P2 and B. 
Therefore it can be neglected in our calculations. 

Thus far we have assumed the factorization approximation to evaluate the decay ampli- 
tudes. It is known in the QCD factorization approach [55] that factorization is justified 



in the heavy quark limit where power corrections of order 1/mg and l/mj can be ne- 
glected. Beyond the heavy quark limit, factorization is violated by power corrections 
which in general cannot be systematically explored. Nevertheless, some of them are 
calculable. For example, in the B decays we have included the terms proportional to a% 
and as which are of order A /m& but chirally enhanced. Final-state interactions which 
have been neglected so far are also of order K/rriQ. The decay — ► (7r + 7r + 7r~)NR pro- 
ceeds only through the iy-annihilation process. Even if the short- distance contribution 
to the weak annihilation vanishes, it may receive sizable long-distance contributions 
via final-state rescattering. The preliminary FOCUS measurement of this mode may 
provide the first indication of the importance of final-state interactions for the weak 
annihilation process in nonresonant D decays. A precise measurement of this mode can 
test the validity of the applying the factorization picture to the nonresonant three-body 
decays. 
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